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ABSTRACT 

A diffusion-limited aggregation process, in which clusters coalesce by means of 
3-particle reaction, A + A + A — > A, is investigated. In one dimension we 
give a heuristic argument that predicts logarithmic corrections to the mean- 
field asymptotic behavior for the concentration of clusters of mass m at time 
t, C m (t) ~ m _1 / 2 (log(t)/t) 3/4 , for 1 < m < y/t/log(t). The total concentration 
of clusters, C(t), decays as C(t) ~ ^/log(t)/£ at t — ► oo. We also investigate the 
problem with a localized steady source of monomers and find that the steady- 
state concentration C(r) scales as r _1 (log(r)) 1 / 2 , r _1 , and r -1 (log(r)) -1 / 2 , re- 
spectively, for the spatial dimension d equal to 1, 2, and 3. The total number of 
clusters, N(t), grows with time as (log(t)) 3 / 2 , t 1 / 2 , and t(log(t)) -1 / 2 for d = 1, 
2, and 3. Furthermore, in three dimensions we obtain an asymptotic solution for 
the steady state cluster-mass distribution: C m (r) ~ r _1 (log(r)) _1 $(z), with the 
scaling function <&(z) = z -1 / 2 exp(— z) and the scaling variable z ~ m/ v/log(r). 

P. A. C. S. Numbers: 82.20.-w, 05.40,+j., 02.50.-r, 82.70.-y. 
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I. Introduction 

Diffusion-limited aggregation processes have attracted considerable recent interest in 
many fields of science and technology [1] . Typically, aggregation processes can be described 
by the binary reaction scheme 

. . K(i,j) 

Here Ai denotes a cluster consisting of i monomers, an i-mer, and K(i,j) is the rate at 
which the reaction between an i-mer and a j-mer proceeds. Much of the understanding 
of the kinetics of binary aggregation processes is based on the analysis of rate equations 
and their exact and scaling solutions [1]. For sufficiently low dimensions, the diffusion 
mechanism is not efficient enough and fluctuations in the densities of diffusing reactants 
result in dimension-dependent kinetic behavior at long times [2]. 

In a view of the richness of the kinetic behavior observed in the bimolecular model, it 
is of interest to investigate more complicated many-particle diffusion-limited aggregation 
processes. In the present study, we focus on the 3-particle reaction scheme, 

A, + Aj + A k K( ^ k) A i+j+k , (1) 

both for homogeneous and inhomogeneous situations, and outline a generalization for the 
n-particle case. 

We will study the simplest n-particle aggregation process for which both the reaction 
rates and diffusion coefficients do not depend on masses of clusters, K(i\, ■ ■ ■ , i n ) = const 
and Dk = const. A notable feature of this process is that it reduces to a simple chemical 
reaction scheme, nA — > A, if one considers only the concentration of clusters. Hence in 
the rate equation description, the concentration, C(t), obeys 

f = -^"< < 2 > 

with 7 being the rate constant. In the long-time limit, the concentration behaves as 

C^Hn-l)*]- 1 ^. (3) 

However, the mean-field rate equation approach provides an accurate description of the 
kinetics only above the upper critical dimension, d c ; when d < d c , the mean-field theory 
does not predict correct long-time behavior. For the aggregation model with constant 
reactivities and diffusivities the upper critical dimension is known: d c = 2/(n — 1) [3- 
5]. Thus, for the binary reaction in one dimension the kinetics is anomalous and the 
concentration decays as t -1 / 2 (see, e.g., [2] and references therein) while the mean-field 
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result is C(t) ~ £ _1 . On the other hand, for n > 4 the mean-field answer (3) gives a correct 
asymptotic description of the kinetics in one dimension. The 3-particle case is marginal in 
one dimension and hence a logarithmic correction has been expected [3-5] . After a number 
of attempts [3-6] , the logarithmic correction of the form C (t) ~ ^log^/t has indeed been 
observed in a very recent study [7]. 

In the next Section II, we explore the diffusion-limited 3-particle aggregation process 
(1) in one dimension. We justify heuristically the appearance of the logarithmic correction 
for the concentration of clusters C(t). Moreover, we obtain a complete scaling description 
of the cluster-mass distribution function: C m (t) ~ (log(t)/t)$(z), with the scaling function 
<&(z) = z -1 / 2 exp(— z) and the scaling variable z ~ m\J\og(t) /t. 

In the last Section III, we examine the 3-particle aggregation process with a spatially 
localized source of monomers. We show that the system reaches the steady state for 
arbitrary spatial dimension d. We present evidence of two critical dimensions, d c = 1 as in 
the homogeneous system and d c = 3 which demarcates the pure diffusive regime (d > 3) 
from the diffusion-reaction regime (d < 3). In three dimensions, we derive a complete 
asymptotic solution for the steady state cluster-mass distribution. 

II. Aggregation process on a line 

Consider a linear lattice on which point clusters undergo a random walk. If n clusters 
happen to occupy the same lattice point, they aggregate irreversibly into a single cluster 
whose mass is equal to the sum of the masses of n parent clusters. This is the n-body 
particle coalescence model (PCM) [8]. The binary PCM in one dimension has been solved 
by Spouge [9], see also [10-14] for other exact solutions of several generalization of the 
PCM. Note that we must keep the lattice spacing Ax finite, even in one dimension, since 
otherwise the reaction will be absent in the n-particle model for all n > 3. The possibility 
of passing to the continuum limit in the binary PCM significantly simplifies the analysis 
(see, e.g., Refs. [10,11,14]), while for n > 3 the spacing Ax, which may be considered as 
the size of particles, appears in the final results. 

A simple heuristic argument explains the logarithmic corrections for the 3-body PCM. 
Let T be a typical time between successive 3-particle collisions in which cluster takes place. 
So the reaction rate is proportional to C/T, 

dC C , A . 

* * ~T (4) 

We will estimate T as follows: Consider a reference frame at rest with respect to an 
arbitrary "target" particle. When two other particles will be at the origin simultaneously, 
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the target particle will die. Now consider all possible pairs of original particles. In the 
following, we will refer to these pairs as imaginary particles. For any pair, let us choose the 
location of one partner as the x-coordinate and the other location as the y-coordinate of the 
corresponding imaginary particle. Thus, we map the original diffusion process on the ID 
lattice onto a diffusion-like process on the 2D square lattice with the same lattice spacing 
Ax. Although imaginary particles do not undergo a simple random walk, we shall assume 
that the asymptotic behavior of this diffusion process is similar to the one encountered in 
the simplest 2D random walk. 

Now let us estimate the collision time T by considering an idealized 2D simple random 
walk of imaginary particles. The target particle will die when some imaginary particle will 
arrive at the origin. The density P(r, t) of imaginary particles is governed by 

^>=MP M ), (5) 

with the initial condition 

P(r,t = 0) = C 2 (6) 

indicating the obvious fact that the density of imaginary particles is just the square of the 
density of original ones, and with the adsorbing boundary condition 

P(r = Ax,t) = 0, (7) 

A simple way to find an approximate solution of Eqs. (5)-(7) is to use a quasistatic 
approximation (see, e.g., [15]). In this approximation one solves the steady diffusion equa- 
tion and accounts for the time dependence by a moving boundary condition. This very 
simple approach often gives asymptotically exact results (see, e.g., [16]). In the present 
problem, the quasistatic approximation yields 

P(r,t) = C 2 l0g( ^ Ax) . (8) 

Now the collision time T may be evaluated by computing the flux to the origin and 
then by equating the flux to the unity: 

J dP(r = Ax,t), 

J 2nD ^ > dt = 1. (9) 
o 

This gives the final estimate, T ~ log(l/C Ax y/2n)/2n DC 2 . Substituting this result into 
(4), we arrive at the closed- form approximate differential equation for the concentration 
of clusters: 

dC 2-kDC . , 

-r- (10) 

dt \og{CAxV^) 
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Solving this equation at the long time limit one gets 



This result confirms previous suggestions [3-6] of possible logarithmic corrections to 
the mean-field power-law decay. Moreover, it is in a complete agreement with recent 
simulational results [7]. A similar logarithmic correction has been found in the binary 
PCM in two dimensions, C(t) ~ log(Dt/Ax 2 )/Dt [17,8,18]. From these findings for 2- and 
3-body PCM in their critical dimensions, d c = 2 and d c = 1, respectively, one can expect 
the appearance of similar logarithmic corrections for the n-body PCM at the corresponding 
critical dimension d c = 2/(n - 1): C(t) ~ (log(Dt / Ax 2 ) / Dt) 1 ^ 71-1 ^. Very recently, a 
similar result has been established by a renormalization group calculation for the n-particle 
annihilation process, nA — > 0, see [19]. Notice also that for d < d c dimensional analysis 
suggests a power-law decay with the exponent dependent on d and independent on n, 
C(t) ~ (Dt)- d / 2 . 

Let us now consider the 3-body aggregation process (1). On the mean-field level, the 
process is described by the (generalized) Smoluchowski equation. For constant reactivities, 
the Smoluchowski equation has the form: 

^ = CiCjCk — 3C m C 2 , (12) 

i+j+k=m 

with C = Notice that in Eq. (12) we set the reaction rate equal to 1 by an 

appropriate choice of units. 

For the monodisperse initial conditions, Ci(0) = No5n, one can find the solution to 
Eq. (12), see [4]. The densities of even-mass clusters vanish due to our choice of initial 
data, while the odd-mass densities are 

c -c r ( m + 1 / 2 ) U c \ m (U ) 

Cl+2m - Cl T(l/2)T(m + l) V~W ) ' (13) 

where V is the gamma function. For example, the concentrations of monomers and clusters 
are given by 

d = N (l + AN 2 t)~ 3 /\ C = N (l + AN 2 t)- 1 ' 2 . (14) 

Note that in the long-time limit the cluster-mass distribution (13) exhibits a scaling 
behavior of the form 

iVo y^rtz 
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with the scaling variable z, z = mC/N . Equations (13) and (15) show that the con- 
centration decays as Ci +2m ~ m~ 

-l/2 t -3/4 for 1 

<C m <C \ft. The exponents describing 
the scaling behavior might be obtained without appealing to the complete solution. For 
example, the exponent describing the time dependence, Ci + 2 m ~ t -3 / 4 , can be derived by 
solving Eq. (12) directly for small m. 

In analogy with the analysis of the chemical reaction scheme, 3A — > A, at the upper 
critical dimension one can suppose that the aggregation reaction process (l)atd = d c = l 
should be described by improved rate equations: 

L ^T= E CiCiCk-SCmC?, (16) 

i+j+k=m 

with L being the logarithmic factor, 

L = log(l/CAxV2^)/(nD). 

Upon summing equations (16) over all m, one can reproduce Eq. (10) for the cluster 
concentration, thus providing a useful check of self-consistency. 

By applying a generating function technique, one can find an asymptotic solution to 
the modified Smoluchowski equation (16). The concentrations are given formally by the 
same expression (13) but with the modified monomer density 



d ~ N ~ 1/2 



log(2Dt/Ax- 
8nDt 



3/4 



(17) 



Both the scaling form and the scaling variable are identical to the corresponding mean- 
field result (15), although the normalization factors and widths of distributions differ by 
logarithmic factors. 

A straightforward generalization to the n-body PCM shows that at the critical di- 
mension, d = d c = 2/(n — 1), the cluster-mass distribution exhibits the scaling behavior of 
the mean- field form [4]: 



C 2 



n-2 



C 1+(n _ 1)m (t) — 1 z «-ie 2 , (18) 



with the scaling variable z, 

z = mC/N Q ~ m(log(L>t/Ax 2 )/L>t) 1/(n - 1) . 
In particular, the concentration of momomers decays as 



Ci(t) 



log(Dt/Ax 2 ) 
Dt 



n/(n-l) 5 



(19) 
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Thus for n-body PCM we have obtained the scaling description of the cluster-mass 
distribution function. Our approach is heavily based on the assumption of equal diffusiv- 
ities of the reactants. [The second assumption of equal reactivities directly follows from 
the assumption of equal diffusivities since "sizes" of clusters in the PCM are the same]. 
However, from our results one can forecast some features of the general case. Consider, for 
example, the 3-body PCM with the cluster diffusion coefficients D m related to their mass 
m by D m ~ m~ 5 . Since the typical mass of the cluster grows inversely proportional to 
the total concentration of clusters, m tyv ~ C _1 , one can insert an estimate for the typical 
diffusion coefficient, D typ ~ m typ ~ C S , into expression (11) for the total concentration of 
clusters. This yields the following heuristic estimate for C(t), 



C(t) 



t 



1/(2+5) 



(20) 



but such an approach cannot predict the scaling form of the cluster-mass distribution 
function. 

III. Aggregation process with a localized steady source 

In this section we consider n-body PCM with a steady spatially localized monomer 
input. We again focus on the 3-particle case and also write some final results for arbitrary 
n. Let the source of monomers is placed at the origin of rf-dimensional space and J is the 
strength of the source. Within the continuum approximation, the concentration of clusters 
at time t and at distance r from the source, C(r, t), satisfies the reaction-diffusion equation 

^p^- = DAC(r,t)-DR+JS(r), (21) 

Here DR is the reaction term: R = C 3 , -C 3 / log(C), andC 1+2 / d for d > l,d = 1, and 
d < 1, respectively. For d > 1, the reaction term has the mean-field form; for the marginal 
case d = 1, the form of the reaction term has been derived in Section 2 (we ignore numerical 
factors and set Ax = 1 by an appropriate choice of the length scale); and for d < 1, one 
can choose the reaction term of the form C 1+2 / d since it gives the correct long-time decay 
for the homogeneous system, C(t) ~ (Dt)~ d / 2 . 

For the binary PCM, a qualitative investigation of similar reaction-diffusion equation 
has been performed in Ref. [12]. Following the same line of reasoning, let us assume that 
the system reaches the steady state and try a solution of the power-law form. By inserting 
such a form to the governing equation one finds that dimension d = 3 also plays a role of 
critical dimension in the present problem. For d > 3, clusters do not interact far away from 
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the source and the concentration decays as r -( d - 2 ), i.e., as in the limit of no reaction. The 
reaction leads to the renormalization of the strength of the source but does not change the 
behavior qualitatively. In other dimensions, the reaction is relevant at all scales, and far 
away from the source C(r) behaves as 



C(r) 



The critical cases d = 3 and d 



r-^logCr))- 1 / 2 , 
r-^logtr)) 1 / 2 , 



if d = 3; 
if 1< d < 3; 
if d = 1; 
if d< 1. 



(22) 



1 have been treated separately. Since we expected 
logarithmic corrections to the power-law behavior C(r) ~ r _1 for 1 < d < 3, we tried 
a solution of the form C(r) ~ r _1 (log(r)) !y , inserted this form into Eq. (21), and finally 
found the exponent v by asymptotically equating the most significant terms. 

Notice that the system with a localized source but without reaction reaches the steady 
state only for sufficiently large spatial dimension, d > 2. On the other hand, the system 
with reaction reaches the steady state for any d. Notice also the appearence of logarithmic 
corrections to the power-law behavior in two critical dimensions, d c = 1 and d c = 3: the 
former results from the logarithmic factor in the reaction term, while the latter reflects the 
fact that at d = d c , the reaction just becomes relevant. 

Since clusters perform a random walk and since the source was turned on at t = 0, 
clusters will propagate diffusively up to the distance of the order \ft. Therefore at r < \ft 
the concentration of clusters approaches to the steady state limit given by Eq. (22); on the 
other hand, C(r) — > rapidly for r > y/i. Hence, the total number of clusters, N(t), may 
be estimated from the relation N(t) ~ f^f* r d ~ 1 C{r)dr. This yields 



N(t) 



taog(t))- 1 / 2 , 

f (d-l)/2^ 
(log(t))3/2, 

I log(t), 



if d > 3; 
if d = 3; 
if 1< d < 3; 
if d= 1; 
if d < 1. 



(23) 



A straightforward generalization to the n-body PCM shows that, again, two critical 
dimensions d c and d c demarcate different behaviors: d c = 2/(n — 1) and d c = 2n/(n — 1). 
For the steady-state cluster concentration one finds 



/ r -(d-2) 



C{r) 



; (log(r)) 



-d c /2 



< r 



t ^ c (log(r) )^ c / 2 , 
-d 



if d > d c ; 
if d = d c ; 
if d c < d < d c ; 
if d = d c ; 
if d < d c ; 



(24) 
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while the total number of clusters scales as 



t, if d > d c ; 

t(log(t))" dc / 2 , iid=d c ; 
N(t) ~ I t( d - d ^/ 2 , if d c < d < d c ; (25) 

(l0g(t)) dC /2 ? if d = dc ; 

k log(t), if d < <i c . 

Consider now the behavior of the steady state concentrations, C m (r), for the 3-body 
PCM in the most interesting three dimensional case. The concentrations C m (r) obey the 
reaction-diffusion equations 



{^2 + l^) C m+ E C i C j C k -3C m C 2 = -(J/D)8 ml 8(r). (26) 

i+j+k=m 



From these equations one can subsequently find asymptotic solutions for the total cluster 
concentration 

C(r) ~r- 1 (41og(r))- 1 /2 5 (27) 
for the concentration of monomers 

dCr^r-^logCr))- 3 ' 4 , (28) 

etc. These results suggest the ansatz 

C m (r) = F m {p)/r, (29) 

with p = log(r) for the behavior of C m (r) for general m. 

Substituting this ansatz into the governing equation, we obtain 

dFm - E FiFjF k -3F m F 2 , (30) 



with 



dp dp 2 

i+j+k=m 



F = YT 1 F ^{P) = rC(r) ~ (4P)" 1 / 2 . 



Far away from the source, one can omit the second term in the left-hand side of this 
equation. After replacement of all F m (p) by C m (t), the resulting approximate equation for 
source-induced aggregation at the steady state becomes identical to Eq. (12) for irreversible 
aggregation, with p, p = log(r), playing the role of time t. For the latter problem, the 
cluster- mass distribution C m (t) approaches the scaling form not only in the monodisperse 
case but also for arbitrary (rapidly vanishing) initial conditions [4]. In general case, the 
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scaling solution still has the form (15), where parameter N is equal to the total mass of 
the system, N = X^m=i m Cm(t) [4]- Since for the former problem the density of mass 
M(r), M(r) = Ylm=i m( ^m{f)-, satisfies the Laplace's equation, 

one gets M(r) = J/AixDr. Therefore Xlm=i m -^m(p) = rM{r) = J/AitD will play the role 
of iVo in the scaling solution. Combining all these findings, we finally obtain that in the 
scaling limit 

m — > oo, r — > oo, z = 2itDm/J^\og(r) = finite, (32) 
the steady-state cluster-mass distribution reaches the scaling form: 

^ / \ tt-D exp(— z) 

Ci+2m r )~ Py J . 33 

Jr log(r) s/ttz 

Thus in the steady state in three dimensions, the typical mass of clusters at distance r 
from the source grows unexpectedly slowly, m typ ~ ^/log(r). 

From the steady-state (33), one can find the following limiting behavior 

C 1+2m (r) ~ ^m-VV-HlogW)- 3 / 4 , (34) 

for m J y/\og(r) / D . Notice that in the steady state the total cluster concentration C{r) 
does not depend on the strength of the source, J, while all other concentrations C m (r) do 
depend on J. 

For the general n-particle PCM, a complete asymptotic solution for the steady-state 
cluster-mass distribution may be found at the critical dimension d = d c following the 
procedure used for the 3-body case. By applying the ansatz C m (r) = F m (p)/r dc one 
can recast the steady-state reaction-diffusion equation to the the n-particle Smoluchowski 
equation. Making use of the scaling solution of the latter equation (18), one obtains 

[" "1 — d-c _ n -2 _ 

Ci+(„_i) m (r) ~ ^rlog(r)j z -'e 2 , (35) 

with z ~ mlog(r) -1 /( n-1 ) being the scaling variable. 

In conclusion, we have obtained the complete asymptotic solution for the steady-state 
cluster-mass distribution for the n-particle aggregation model with a spatially localized 
source of monomers in the critical dimension d = d c = 2n/(n — 1); the binary case cor- 
responds to d = 4 while the ternary case corresponds to d = 3. The solution exhibits an 
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unexpectedly slow logarithmic growth of the typical mass of clusters versus the distance r 
from the source: m typ ~ \og(r) 1 ^ n ~ 1 \ For sufficiency large dimensions, d > d c , we have 
found that clusters do not interact far away from the source and hence concentrations de- 
cay diffusively as r -( d ~ 2 ). Therefore for d > d c the reaction leads to the renormalization of 
the strength of the source but does not change the behavior qualitatively. For sufficiently 
small dimensions, d c > d > d c: the steady-state cluster- mass distribution is expected to be 
of the scaling form C m (r) ~ r _/3 $(mr~ Q ). From two known moments of the cluster-mass 
distribution 

C(r) = Y°° C m (r)~ r - d ° 

' 'm=l 

[see Eq. (21)] and 

M(r) = V°° mC m (r) ~ r"^" 2 ), 

one can find the exponents a and (3, a = d c — d and (3 = d c + d c — d. However the scaling 
function $(mr~ Q ) satisfies a rather difficult integro-differential equation which I could not 
solve. Further progress is also possible in one dimension where exact solutions for the total 
density of clusters and for the densities of clusters of small mass may be readily found and 
behaviors even in the vicinity of the source may be investigated. 

I am grateful to D. ben-Avraham for informing me on his results [7] prior to publi- 
cation and to E. Ben-Nairn, F. Leyvraz, and S. Redner for useful discussions. I gratefully 
acknowledge ARO grant #DAAH04-93-G-0021 and NSF grant #DMR-9219845 for partial 
support of this research. 
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